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Abstract 

Questions concerning the Bose-Einstein (BEC) and Fermi-Dirac (FDC) 
correlations of hadrons produced in e~^e~' collisions are discussed. Among 
them the emitter dimension r as a function of ^/s^ and the hadron mass, 
the extension of the BEC by including isospin invariance and the proposed 
relation between r and the inter-atomic separation in Bose condensates. 



1 Introduction 

The present report covers a few of the open questions which are related to the 
Bose-Einstein and Fermi-Dirac Correlations (BEC and FDC) of hadrons in e~^e~ 
annihilations. The topic which are dealt here include the dependence of the 
hadron emitter dimension on the ^/s^ and on the mass of the hadron; the in- 
terpretation of the 2-dimensional BEC analysis results obtained in annihilations 
as compared to those found in heavy ion collisions and finally the status of the 
so called Generalised BEC. The issue whether one may relate some features of 
the BEC at high energies with those present in Bose condensates is also briefly 
discussed. 



2 Dependence of the dimension r on the y^s^ 

As is well known the Bose-Einstein interferometry method rest on the fact that 
the density of identical boson pairs is enhanced when they are emitted near in 
momentum and phase space. Assuming the emitter to be a sphere of a Gaussian 
distribution, the correlation of pairs of identical bosons can be parametrised by 
C2{Q) = 1 -l- X2G~^ ^ where = — (pi — P2Y and pi and p2 are the four 
momenta of the two bosons. An early compilation of the r^^ values measured 
in heavy nuclei (AA) reactions is seen in Fig. 1 to be rather well described by 
r = 1.2 A^^^ fm, the known nucleus radius dependence on the atomic number A. 
This behaviour is taken to be the result of the increase of pion sources as A gets 

^Presented at the ISMD2003, 5-11 September 2003, Krakow, Poland 
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larger. As for the Ttttt values obtained from the e~^e~ annihilation analyses, they 
are seen in Fig. 1 to be essentially independent of ^/s^ as it is generally accepted 
that the small deviations from a constant value are due to the non-uniformity of 
the chosen experimental procedures. Since e.g. the average hadron multiplicity 
increases with energy this fiat r distribution poses a theoretical challenge. 
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Figure 1: Left: r-^^^ versus A^/^ found in heavy ion collisions. The line represents 

versus ^ See found in e~'~e annihilations. The line is the 
expectation of a hadron jet source model (see text). 



In an effort to understand the origin of this flat r distribution we adopt here, for 
the e~^e~ annihilations, the hadron jets as the pion sources and further assume 
that only pion pairs emerging from the same jet may be correlated. Now the effect 
of multi-hadron sources on particles' correlation strength has been studied by 
several authors [1] in the frame work of the factorial cumulants moments Kg where 
g = 2 in the case of pion-pair correlations. Specifically, it has been shown that the 
cumulant is diluted by a factor Dg, so that = Dr,-KP^^ where D„ is defined as 



- . , ^ /y. q""^^ + Ar^^i//e^e"t) . Here N^^""^ and Nf^f^''^''* denote respectively 
the number of q pions combinations coming from the same and from different 
hadron jets. Now for g = 2 there exists the relation C2{Q) = 1 + A2e~'^ = I+K2 
which can further be extended to include the dilution factor D{= D2) so that 
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The dependence of Q in this expression is then removed by integrating the quan- 
tity C2{Q) — 1 over Q from zero to infinity to finally obtain 
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To apply this formula to the e+e^ hadrons reaction there is the need to 
evaluate the dilution factor D as a function of ^/s^. As is well known the de- 
termination of the number of hadron jets in a particle reaction does depend on 
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the chosen jet identification method [2] . Here for simphcity we will assume that 
the e^e~ annihilation proceeds via three hadron jets, two coming from the quark 
and anti-quark pair and the third from a radiated gluon^. The total ^J's^ is then 
divided in the following way. The first quark jet gets on the average the energy 
^/s^/2 and thus the second quark jet has the energy of — Egiuon- The final 

number of outgoing charged pions assigned to each of the three jets is determined 
from the known total average charged hadron multiplicity dependence on ^/s^ 
and the known average charged multiplicity of quarks and gluon jets parametrised 
as a function of their jet energy [2]. From this pion division the dilution factor 
D2 is estimated and the resulting r dependence on ^/s^ is shown in Fig. 1 by the 
continuous line which is normalised at 40 GeV while setting \d — \d=\- As seen, 
the r behaviour extracted from this simplified model, which rises by only ~3.5% 
from 10 to 90 GeV CM energy, already describes nicely experimental situation. 



3 The r dependence on the hadron mass 

The very high statistics of the Z° hadronic decay events accumulated by the LEP 
experiments and the extension of the BEG analysis to baryon pairs via the FDG 
method opened the way to study the emitter dimension function of the 

hadron mass, from pions to A baryons. The results of these analyses are shown 
in Fig. 2 together with two recent, L3 and OPAL, measured r^o^o values which 
arc seen to be inconsistent and therefore are presently disregarded. The emitter 
dimension r{m) is seen to decrease with the hadron mass, from ~0.75 fm at ttItt 
to ~0.15 fm at mA- It has further been shown [4] that this behaviour can be 
derived from the Heisenberg uncertainty relations which yield 

r{rn) — c^/A^h/^/m . (1) 

Choosing At = 10^^"^ second, as a representative time scale of the the strong 
interaction sector, r(m) is shown in Fig. 2 by the continuous line. An almost 
identical expectation for r{m) is also forthcoming from the Local Parton Hadron 
Duality and a general QGD potential [4]. On the other hand, the decrease of r 
with the hadron mass poses a challenge to a large variety of hadronisation mod- 
els, including the Lund one, which expect r to increase with the hadron mass. In 
particular at present, there is no satisfactory explanation for the very small r\ 
value. 

In heavy ion collisions part of the emitter survives after the emission of the iden- 
tical bosons used in the BEC analysis. This is not the case in e~^e~ annihilation 
where he emitter collapses with the emission of its hadrons used for the BEC 
analysis. Consequently one may try and estimate an approximate minimum en- 
ergy density of the emitter of the two hadrons by the relation ^'density = '^ttih/V — 

^Even at the Z° the fraction of events with more than three hadron jets is very small if one 
uses e.g. ycut > 0.01 in the Durham jet identification scheme [3]. 
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Figure 2: Left: The r dependence on the hadron mass in hadronic decays. Right: 
The minimum energy density of the hadron emitter [1] versus the hadron mass. 

6m/i/(47rr'^) which is shown in Fig. 2. Whereas the energy density for the pion 
and kaon emitters still lies within a reasonable range, the energy densities of the 
protons and A's emitters are lying in the vicinity of 100 GeV/fm^, way above the 
values anticipated from the current models for hadron production. 

4 The 2-dimensional BEC analyses 

In recent years the correlations of pions emerging from heavy ion collisions and 
from e~^e~ annihilations have been analysed in terms of the 2-dimensional BEC 
analysis. This analysis is carried out in the Longitudinal Centre of Mass System, 
the description of which can be found e.g. in Ref. [1]. In this system the direc- 
tion of the event thrust is referred to as the longitudinal direction and the two 
other directions are referred to as the outgoing and side axes. To these directions 
one associates three four momentum differences Qz, Qout and Qside which in a 2- 
dimensional analysis are reduced to Qz and the transverse momentum difference 

Qt = \jQ'iut + Qlide- ^his method the emitter is not anymore restricted to a 
spherical configuration but can also assume the shape of an ellipsoid. Further- 
more, one has an additional variable, the so called transverse mass iriT-, defined as 
rriT = 0.5 X [^m? + + ^^TiM-pl^) where the pxi are the transverse momenta 
of the two hadrons. From a fit to the data one is able to determine the longitudi- 
nal tl and the transverse tt dimensions. The dependence of on ijit are shown 
in Fig. 3 for e^e~ annihilations, as measured by DELPHI at the Z° mass, and 
for S+Pb collisions at 200 GeV/A as measured by the NA44 collaboration. First 
to note is that riimT) coming from e^e~ annihilations is very similar to r(m) 
shown in Fig. 2 which was obtained from the 1-dimensional BEC analyses of the 
Z° hadronic decay. Moreover it has indeed been shown [5] that, as in the case 
of r(m), also r2,(mr) can be deduce from the Heisenberg uncertainty relations to 
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yield 

riinT-T) ~ At h/y/rriT (2) 

which is the Eq. (1) expression where m and r(m) are substituted respectively 
by rriT and riimT). The dependence of tl on l/^/rnr is repeated also in the 
heavy ion collisions with the difference that the proportionality factor is equal to 
2.0 as compared to 0.354 found for e~^e~ annihilations. To note is that the ratio 
2.0/0.354 is not far from the ratio of the average radius of the S and Pb nuclei 
to the radius of ~0.85 fm obtained from the 1-dimensional BEC of the hadronic 
Z° decays. At this point the obvious question arises as to what extend one is 
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Figure 3: The r as a function of my. Left: DELPHI preliminary results from e~^e~ — > 
hadrons. The line represents r = 0.354/^mT. Right: NA44 results from S+Pb 
collisions. The line represents r = Ij ^pmr- □ □□ 

permitted to extract from the 2-dimensional BEC analyses information on the 
underlying dynamics of heavy ion collisiofg when r2,(mr^e?iaves essentially the 
same as the in e^e~ annihilations which are void of nuclear matter. 



5 Generalised Bose-Einstein Correlatior^ (GBEC) 

In analogue to the generalised Pauli principle, which extends that principle fronP 
two identical nucleons to the proton-neutron system, one may consider an inclu- 
sions of the isospin invariance to extend the BEC to a generalised BEC (GBEC). 
Such an extension has been proposed by several authors [1] with th£Ikim to es- 
tablish BEC relations between e.g. the tt^tt^ and tt^tt" systems. One relation of 
this kind, which was proposed by [5], states thax 

^jP{%o ^ (7r±7r±)X] = ^P[^o ^ (7r±7rO)e.en I X] 

X X 

Here P stands for probability, Iq stands for an 1=0 initial state, like T ^ bb ^ 
hadrons, and X stands for the rest of the final state associated with the specified 
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two pion system. The GBEC validity at high energy physics is still an open 
question due to the missing of an experimental verification. 

6 Conclusions 

The r(m) behaviour should be investigated also in other than e'^e~ particles 
reactions and hopefully the GBEC applicability to hadron production will soon 
be tested experimentally. In trying to gain an insight into the r(m) behaviour 
and the energy density of the hadron emitter it may be instructive to turn to 
other phenomena related to the Bose-Einstein statistics. One of them are the 
Bose condensates where it was shown [5] that at equal very low temperature 
their inter-atomic separations arc proportional to 1 / yjmatom and not to their 
dimension. Moreover, in as much that the Heisenberg relations are applicable to 
condensates, as most of their atoms are in the energy ground state, the inter- 
atomic separation formula can be transformed to Eq. (1) so that a comparison 
with the BEG dimension is not unreasonable. This comparison coupled to the 
behaviour of r{m) and the high emitter energy densities may point to the need 
to re-examine the interpretation given to the BEG derived dimension r. 
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